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ABSTRACT 

O 

We present Very Large Array observations with 80 milliarc second resolution (~ 9 

pc) of the recently discovered Galactic-analog H 2 masers in the Antennae interacting 
galaxies (NGC 4038/NGC 4039; Arp244). Three regions of water maser emission are 
detected: two in the "interaction region" (IAR) and the third ~ 5.6" (> 600 pc) west 
of the NGC 4039 nucleus. The isotropic H 2 maser luminosities range from 1.3 to 7.7 
L Q . All three maser regions are mostly obscured in the optical/near- infrared contin- 
uum, and are coincident with massive CO-identified molecular clouds. The H 2 maser 
velocities are in excellent agreement with those of the molecular gas. We also present 
archival VLA 3.6 cm data with ~ 0"28 (~ 30 pc) and ~ 0"8 (~ 90 pc) resolution to- 
ward the maser locations. All three maser regions are coincident with compact 3.6 cm 
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radio continuum emission, and two are dominated by thermal ionized gas, suggesting 

the presence of natal super star clusters containing the equivalent of a few thousand O 

stars. We also present detailed comparisons between the radio data and existing HST 

ACS (optical) and NICMOS (near-IR) data and find that both maser regions in the 

IAR are also associated with Paa emission and neither source is detected shortward 

of 2 ^m. These results highlight the potential of using Galactic-analog H 2 masers to 

pinpoint sites of young super star cluster formation with exquisite angular resolution. 

m 

Subject headings: galaxies: interactions — galaxies: ISM — galaxies: star clusters: 
general — galaxies: starburst — masers — radio lines: galaxies 
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1. Introduction 

Water masers are found in the vicinity of ~ 70% of infrared brig ht (100 /mi > 1000 Jy 



and 60 /im > 100 Jy) ultracompact Hn (UCHII) regions in our Galaxy (IChurchwell et alJI 1990c 
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Kurtz & Hofnerll2005|) . providing excellent signposts for active star formation. The term "kilo- 
maser" was coined to describe extragalactic H 2 masers with luminosities comparable to the 
brightest star formation H 2 masers in our Galaxy (e.g. W49N; L H2 o ~ 1 L©). Kilomaser isotropic 
luminosities (L H2 o < 10 L Q ) are much lower than the more widely studied " megamasers" found in 
the nu clear regions of AGN (i.e. NGC 4258), with luminosities up to 10 4 L (IBarvainis & Antonucci 
2005J). Since "kilomasers" can in principle be either of nuclear (i.e. amplifying a background 
AGN) or star formation origin, we will use the term "Galactic-analog" (GA) H 2 maser to refer 
to non-nuclear masers with kilomaser luminosity. Natal super star clusters are thought to contain 
1000s of UCHII regions, therefore GA-H 2 masers may help pinpoint sites of young extragalactic 
cluster formation, though there have been few searches with the required sensitivity. Until re- 
cently GA-H 2 masers had been unambiguously detected toward only six nearby galaxies: LMC, 
M82 , IC 342, IC 10, M33, arid NGC 2146 with i sotropic luminosities in the range 0.005 to 4 
L Q (fWMteoak & Gardnerll 1986c iHenkel et a/j2005t ICastangia et a/.ll2008l and references therein). 
The remaining few known H 2 kilomasers are either associated with AGN activity or their origin 
is ambiguous. 

We recently conducted a Greenbank Telescopaj (GBT) search for H 2 masers toward four 
nearby starburst galaxies (3 < D < 22 Mpc), known to harbor natal super star clusters and without 
known AGN (or at worst low luminosity AGN) down to a sensitivity level sufficient to detect strong 
GA-H 2 masers. Positive detect ions were found for all four galaxies (the Antennae galaxies, He2- 
10, NGC 4214, and NGC 5253: barling et al.l2008( ). suggesting that GA-H 2 maser emission may 
be common in starburst galaxies. The unusual success of this mini-survey arises from its order of 
magnitude greater sensitivity compared to most previous single dish surveys for either kilomasers 
or megamasers. 

The strongest masers were detected toward the Antennae (D~ 22 Mpd_|) interacting galaxies 
(NGC 4038/NGC 4039; Arp244). Indeed, this violent merger is an ideal laboratory to study how 
H 2 maser emission behaves in regions of extreme star formation. In this paper we report sensitive 
high angular resolution Very Large Array (VLA) water maser observations of the Antennae galax- 
ies, along with archival 3.6 and 6 cm continuum, CO, and optical/near-IR data in order to pinpoint 
the maser locations and determine the nature of the emission. 



'The National Radio Astronomy Observatory operates the GBT and VLA and is a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universities, Inc. 



2 Antennae distance estimates range from 13.8 ± 1.7 Mpc dSaviane et al\ 120081) to 22 ±3 Mpc dSchweizer et al. 



2008); we adopt 22 Mpc. 
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2. VLA Observations and Data Reduction 

We observed the (616 - 523) ortho-water maser line at 22.23508 GHz toward the Antennae 
galaxies in the VLA A-configuration using fast- switching and reference pointing. The full width 
at half power (FWHP) at the Doppler shifted frequency of ~ 22.11 GHz is ~ 2', encompassing 
both nuclei (NGC 4038 and NGC 4039), as well as the interaction region (IAR; see Figure [QO; 
in comparison the GBT primary beam was only 30". In order to cover the velocity extent of 
maser emission detected with the GBT with similar velocity resolution (3.7 km s" 1 ) required three 
correlator settings, each with 2 IFs. We also retrieved A and B-configuration 3.6 cm continuum 
data, as well as BnA-configuration 6 cm data fr om the VLA archive ( the B-configuration data 



were included in the multi-configuration study of LNeff & Ulvest ad 2000). Further details of the 
VLA data are given in Table [Q Data calibration followed standard high-frequency procedures 
in AIPS, including using a model for the brightness distribution of the absolute flux calibrator 
3C 286; absolute flux calibration is good to ~ 5%. The line-free channels from the ~ 22.1 1 GHz 
data were used to estimate the 1.3 cm continuum emission in the UV-plane. The maser line data 
were subsequently imaged in CASA (to facilitate regridding of the six IFs), and primary beam 
correction was applied. The continuum data were imaged in AIPS, and primary beam corrected. 
The astrometric accuracy is better than ~ 0"05, while the relative position uncertainty of the maser 
data is an order of magnitude better). All velocities are presented in the barycenter frame, optical 
definition. 



3. Results 

3.1. H 2 Maser Emission 

As shown on the three-color Spitzer composite in Figure \T^l, three distinct regions of H 2 
maser emission are resolved by the VLA; two are located in the 24 /im-bright IAR between the 
two nuclei, and one is located 5 '.'6 (> 600 pc) to the west of the NGC 4039 nucleus. For ease of 
referencing, these three maser regions are denoted H 2 0-East, H 2 0-SE, and H 2 0-West. We call 
these "maser regions" because at the current resolution ~ 80 mas (~ 9 pc), the observed emission 
is most likely the sum of many individual maser spots. Spectra of the maser emission are shown 
in Figures [Q),c, and d, and the observed properties are given in Table |2] At the current angular 
resolution, no variation of maser position with velocity is detected. The range of maser velocities 
detected by the VLA is in excellent agreement with that of the GBT data, despite the 16 x larger 
area covered in the VLA data. Differences in intensity between the VLA and GBT data can be 
explained by the maser offsets from the GBT pointing center (see Fig.Qi). Using the observed line 
properties (Table© we find isotropic H 2 maser luminosities ranging from 1.3 to 7.7 L , which 



Table 1. VLA Observing Parameters 



Parameter Value 



H2O maser and 1.3 cm (22.1 1 GHz) continuum Observations 

Project (config.) AB 1 304 (A) 

Observing Dates 10, 11, 15 Nov 2008 

Bandwidth 6 X 6.25 MHz 

Velocity Bandwidth 1360 to 1740 km s -1 

Velocity resolution 3.7 km s -1 

Spec, line resolution 108 mas X 61 mas (P.A.= +1.8°) 

Spec, line rms noise 1 mJy beam -1 

Cont. resolution" 350 mas X 220 mas (PA.= 10°); < 1000 kA 

Cont. rms noise 0.14 mJy beam -1 

Archival 3.6 cm (8.46 GHz) Continuum Observations 

Project (config.) AP478 (A), AS796 (A), AA301 (A), AN079 (B) 

Observing Dates 01 Nov, 21 Dec 2004; 09 Mar 2006; 07 Sep 1998 

Bandwidth 2x50 MHz 

Cont. resolution (A) 350 masx220 mas (P.A.= 10°); > 30 kA 

Cont. rms noise (A) 0.032 mJy beam 

Cont. resolution (B) 1"03 X 0"67( P.A.= -2.4°) 

Cont. rms noise (B) 0.013 mJy beam -1 

Archival 6 cm (4.85 GHz) Continuum Observations 

Project (config.) AN074 (BnA) 

Observing date 31 Jan 1997 

Bandwidth 2 X 25 MHz 

Cont. resolution 1 "03 X 0" 67 (PA.= 82. 1 ° ); > 65 kA 

Cont. rms noise 0.022 mJy beam -1 



Note. — A Briggs weighting of robust=0 was used unless otherwise speci- 
fied. UV weighting adjustments are indicated after the resolution; < indicates 
a Gaussian taper and > indicates a cut. 

"Natural weighting was used. 



is on the high side, but not dissimilar to the other galaxies with known GA-H 2 emission (see 
§1). The H 2 luminosity sensitivity limit is ~ 0.6 L (assumes V peak =5 mJy (5a) and /SV FWHM =\Q 
km s" 1 ), brighter than the majority of Galactic water masers. 



3.2. Radio Continuum Emission 

We have created matched resolution ~ 0"8 (~ 90 pc) 3.6 and 6 cm continuum images with 
high signal-to-noise, as well as higher resolution ~ 0''28 (~ 30 pc) 1.3 and 3.6 cm images, albeit 
with significantly less sensitivity. Care was used to match the UV-coverage of the continuum pairs 
at both resolutions. These images are the highest resolution radio images avai lable to date tow ards 
the Antennae (see for example the ~ 1.1" resolution 3.6 cm VLA data of lNeff & Ufvestadl2000l) . A 



complete analysis of these data throughout the Antennae will be presented in a forthcoming paper, 
for now we concentrate on the maser locations. Radio continuum sources associated w ith H 2 Q-East 



and H 2 Q-West are detected in the ~ 0''8 resolution 3.6 and 6 cm images (also see lNeff & Ulvestad 



20001 their sources 4-4 and 1-3 respectively). We have also resolved compact continuum emission 



towards H 2 Q-SE in the ~ Of 28 1.3 and 3.6 cm images (this source is not resolved from the bright 



Neff & Ulvestad! I2000L source 2-1 at poorer resolutions). The radio source associated with H 2 0- 
East is also detected in the higher resolution (poorer sensitivity) 1.3 and 3.6 cm continuum images, 
though H 2 0-West is not. The 3.6 cm continuum contours are shown in Figures [2^,b,c,d (see §4.1). 

The radio continuum sources are denoted CM-H 2 0-East, CM-H 2 0-SE, and CM-H 2 0-West, 
and their observed flux densities and spectral indices (S u o c v a ) are listed in Tabled The spectral 



index for CM-H 2 0-East is flat and in good agreement with lNeff & Ulves tad (2000, their source 4- 



4, a = -0.16 ±0.19). We also find aflat (or possibly inverted) spectrum for CM-H 2 0-SE (detected 
here f or the first time). I n con trast, we find that CM-H 2 0-West has a steep spectrum (a = —0.6^ 4), 



while iNeff & Ulvestad! (12000) inferred an inverted spectrum (their source 1-3, a = +0.38 ±0.19). 



Table 2. VLA Water Maser Properties 



Maser 
Region 


Position 
a (J2000) 8 (J2000) 


V range 

(kms- 1 ) 


'peak 

(kins" 1 ) 


dpeak 

(mJy) 


fSdv 
(mJy km s ' ) 


(£©) 


H 2 0-East 

H 2 0-SE 

H 2 0-West 


12:01:55.4590(2) 
12:01:54.9959 (4) 
12:01:53.1257(1) 


-18:52:45.653 (5) 
-18:53:05.543 (8) 
-18:53:09.805 (2) 


1502 - 1528 
1439 - 1539 
1639 - 1689 


1520.4 
1515.1 
1670.5 


9.3 (0.9) 
7.2(1.0) 
44.3(0.1) 


116(12) 
370 (25) 
690 (25) 


1.3(0.1) 
4.1 (0.4) 
7.7 (0.5) 



Note. — Positions and peak flux densities measured from 2-D Gaussian fits. Isotropic line luminosities were compute d from 
Lfap = (2.3 X 10~ 5 Lq) X D 2 X J Sdv, where the distance D was assumed to be 22 Mpc and f Sdv is in mJy km s~' IHenkel et al\ 

2005). 
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Fig. 1. — (a) Three color Spitzer image of the Antennae galaxies with RGB mapped to 24, 8, & 
3.6 jum. Green + symbols show the locations of the three H 2 maser regions detected by the VLA. 
In panels (b), (c), and (d) the VL A H?0 maser spectr a are shown by the solid histogram, and ~ A" 
resolution CO(1-0) spectra from Wilson et alJ (|2000t ) are superposed as dot-dash lines. 
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This inconsistency can be explained by the fact that these authors use a fitted size for this source 
at 3.6 cm, but only a peak flux density at 6 cm due to a poor Gaussian fit, thus underestimating 
the 6 cm flux density. Since both CM-H 2 0-East and CM-H 2 0-SE have flat radio spectra, and are 
coincident with Pact emission (see§4.1), both are likely to be dominated by thermal Hll regions. 
In contrast, CM-H 2 0-West has a steep spectrum and thus appears to be dominated by non-thermal 
emission; broader radio wavelength coverage is required to disentangle a thermal component for 
this maser region. 

For the two thermal sources, CM-H 2 0-East and CM -H 2 Q-SE , the p roduction rate of ionizing 
photons (<2l> c ) can be e stimated using Equation 2 from Condon! (| 19921) assuming thermal, opti- 
cally thin emission (see I Johnson et al.1 12009). Although the 1.3 cm data would generally be pre- 
ferred (to exclude non-thermal contamination, and to have a higher likelihood of being optically 
thin), we have used the 3.6 cm data to determine Q L>C due to its superior se nsitivity (see Ta ble [p. 
An 07.5V star p roduces an ionizing flux of Q Ly c =10 49 s" 1 (hereafter O*: iLeitherer et al.lll999i : 
Vacca et al.lll996|) . suggesting that the ionized gas associated with the two IAR maser regions is 
equivalent to ~ 2000-5000 O* stars. Note that these Ql } c values could be underestimates if a 
signi ficant fraction of the ionizing flux is absorbed by dust or able to esc ape from the H II region 



(see Joh nson & Kobulnickyll2003 



Reines et al. 2008 



Johnson et ai]l2009h . 



4. Discussion 



4.1. Multi wavelength Comparison of the Maser Environments 

To further assess the nature of the sources associated with the maser emission, we have ob- 
tained existing Hubble Space Telescope (HST) ACS and NICMOS data, along with Owens Valley 
Radi o Observatory (OVRO) CO(l-O), and Submillimeter Array ( SMA) CO(3-2) molecular line 
data dwhitmore et all 1 9991 : Iwilson et all200ol : IPetitpas et alJ2007h . To use these data we must first 
ensure the registration of the different data sets. The astrometry of the VLA data is good to 0"05, 
far superior to the other available data. Because of their shorter wavelengths, the OVRO and SMA 
CO positions are likely less accurate (~ O'.'l). The HST data hav e an astro metric uncertainty of 
~ I". For example, by comparing with the lNeff & Ulvestadl (|2000 ) radio data. IWhitmore & Zhang 
(120021) found a V.'2 offset for their WFPC2 images. For the current analysis we first matched the 
continuum-subtracted Pact emission in the IAR to the 3.6 cm continuum since they have the same 
physical origin. We then applied that astrometric solution to the non-continuum-subtracted Paa 
image, which has many features in common with the other optical/near-IR images, allowing us to 
register and correct all the other HST images. The resulting shifts for the ACS and NICMOS data, 
are -0"97, +1 "58; and +0"49, +0"20, respectively (we also find good agreement with the WFPC2 



Table 3. Radio Continuum Properties 



Source F 6 . 0cm a F 3 . 6cm a F 13cm a a Q LyC c 

(mJy) (mJy) (mJy) (xlO^s" 1 ) 

CM-H 2 0-East d 1.1(0.1) 1.00(0.08) ••• -0.17^i 5.1 

CM-H 2 0-SE e ••• 0.5(0.1) 0.6(0.4) +0.2!^ 1.9 

CM-H 2 0-West d 0.9(0.14) 0.64(0.09) ••• -0.6% ? 

a Flux densities measured from polygon regions; uncertainties are (no. independent beams) 5 x 
3<t plus a 5% calibration uncertainty (added in quadrature); blanks mean source was resolved out 
or too confused. 

b Uncertainties are 68% confidence levels (la); for H 2 0-SE 0.4 mJy is a lower limit for the 
1.3 cm flux density, so this defines the lower bound on its a. 

c Inf erred from 3.6 cm flux densities. 

d Flux densities measured from the ~ 0.8" data. 

e Flux densities measured from the ~ 0.28" data. 
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Fig. 2.— Three color images of (a) the IAR, (b) NGC 4039 nucleus, (c) H 2 0-East, and (d) H 2 0-SE. 
For (a) and (b) RGB is mapped to HST ACS optical filters Ha, 8 14W, & 550M. For (c) and (d) RGB 
is mapped to HST NICMOS near-IR filters Paa, F160, & Fl 10. In all four panels blue + symbols 
show the H?0 mase r locations and the yellow contours show SMA CO(3-2) integrated intensity 
(IPetitpas et al J 120070 . Green contours show VLA 3.6 cm 0''8 resolution continuum emission in 
(a), (b), and (c) with contour levels of 0.05, 0.1, 0.3, 0.6, and 1.2 mJy beam" 1 . Green contours 
show VLA 3.6 cm 0''28 resolution continuum emiss ion in (d) with cont our levels of 0.1, 0.2, 0.3, 
0.6, and 0.9 mJy bea m" 1 . The SGMCs identi fied by lWilson et alj (|200Q[) and the cm- A continuum 
sources identified by iNeff & Ulvestadl (I2000|) are indicated. 
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offset determined by iWhitmore & Zhang||2002|) . We estimate that the final relative astrometry is 
good to ~ Of 2. 

Figures^, and b show three-color optical wavelength HST images of the IAR, and NGC 4039 
nuclear region, respectively, and are superposed with the locations of t h e mas ers, 3.6 cm contours, 
and SMA CO(3-2) integrated intensity contours from IPetitpas et al.1 (|2007l) . Figures |2fc, and d 
further zoom into the IAR to show close-up near-IR NICMOS views of H 2 0-East and H 2 0-SE, 
respectively. Addition ally, superposed on the H 2 spectra in Figures [Tp,c, and d are scaled OVRO 
CO(l-O) spectra from lWilson et al.l (|2000l these data have significantly higher S/N than the CO(3- 
2) data). The velocities of all three maser regions are encompassed by the velocity range of the 
CO(l-O) emission, suggesting that the masers are located within the thermal molecular gas. In- 
deed, H 2 0-East and H 7 Q-SE are coincid ent with super giant molecular complexes (SGMC) 1 and 
4-5, respectively (seelWilson et al.ll2000h . and all three maser regions are coincident with CO(3-2) 
clouds (also see lPetitpas et al.ll2007l) . Additionally, all three maser regions are notably dim at opti- 
cal wavelengths, suggesting high obscuration (see Fig. [2^, b). Below we discuss the properties of 
each region in detail. 



4.1.1, H 2 -East Region 

As shown in Fig. [2k and discussed in §3.2, the H 2 0-East maser is coincident with a massive 
thermal H II region, but is located in an optically uninteresting region between two optically bright 
star forming regions ~ 4" (> 430 pc) to the NE and SW. In the following, we will denote the H II 
region itself as H 2 0-East. H 2 0-East is also coincident with a source of compact Pa« emission 
detected by NICMOS (Fig. |2fc), though it is not detected in either of the NICMOS F110W, or 
F160W continuum images (it is weakly detec t ed in the F220M image, not shown). The H 2 0- 
East cluster is not included in the iBrandl et al.l (120051) WIRC J and Ks band study of the reddest 
Antennae clusters because it is so red that it was not detected at J-band. N otably H 2 Q-East is only 



1"5 west of the 2nd reddest Antennae cluster (ID W i RC =190, A v ~ 10; brandl et all [2005). In 



contrast, the optically bright massive clusters located ~ 4" to the NE and SW of H 2 0-East show 
relatively little extinction (A v =1.3 for the SW sources IWhitmore & Zhang 2002). Unfortunately, 
the resolution of existing mid-IR (A > 2.2 fim) imaging has been too poor to resolve the H 2 0-East 
region from the nearby clusters. 



H 2 0-East has also been detected in Br7 by iGilbert & Graham! (120071 their source "D2") 
using NIRSPEC on Keckll. These authors find a peak Br7 velocity of 1487 km s" 1 and a AV FWH m= 
55 km s _1 , which is somewhat blue-ward of the H 2 0-East maser and CO(1-0) peak velocity of 1520 
km s" 1 ( see Figffp), sugges ting the Br7 originates in outflowing ionized gas in front of the maser 
region. IGilbert & Graham! (|2007l) find that the bright cluster to the SW (their source "D", and 
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Neff & Ulvestadl 120001 source CM-4-2) has a Br7 velocity offset by 100 km s" 1 from H 2 0-East, 
suggesting that these clusters may be separated significantly along the line of sight in addition to 
the ~ 4" (> 430 pc) plane-of-sky distance. 



4.1.2. H 2 0-SE Region 

The H 2 0-SE region is located only 0"5 (~ 60 pc) to the NE of the most massive (~5x 
10 6 M ) and brightest mid-IR to cm-A star cluster ( = WS80, Bl, CM-2-1, WTRC157, Peak 1) in 



the whole of the Antennae ga l axies (see for example IWhitmore & Zhangll2002l : iGilbert & Graham 
2Q07J : lNeff& Ulvestadl bood : brandl et alJ l2005l J2009L we denote this cluster as WS80 for ease 
of referral). Though the optical HST data have more than sufficient resolution, H 2 0-SE is not 
detected. The majority of longer wavelength data do not have sufficient resolution, so that only the 
integrated properties of H 2 0-SE plus the WS80 cluster are availa ble. Indeed, the H 2 Q -SE region 
has been resolved and detected in only one other study - that by ISnijders et al.l (I2006|) using the 
VLT VISIR instrument at 12.81 /mi ([NEII] + continuum, their source "lb"). The H 2 0-SE source 
contributes about 1/5 of the total flux of the combined H 2 0-SE + WS80 regions at this wavelength, 
similar to what is seen in the 0''28 resolution 3.6 cm data. 



Wilson et al.1 (I200QQ find this region to be coincident with the confluence o f two SGMCs, and 



as seen in Fig.QJ, both the maser and CO(1-0) show particularly broad emission. IGilbert & Graham 
J2007I) find a Br7 velocity for WS80 (their source "Bl"), of v peak = 1476 km s" 1 with a FWHM= 
70 km s" 1 . These values are in excellent agreement with the H 2 0-SE maser kinematics (1515.1 
km s" 1 , and 100 km s" 1 ), suggesting that H 2 0-SE and WS80 are in close proximity kinematically 
as well as spatially. 



Estimated ages for t he WS80 cluster ran ge between 1-3.5 Myr (IWhitmore & Zhangll2002l : 
Gilbert & Graham! 120071 ; ISnijders et al.ll2007|) . though it is diffic ult to pinpoint cluster ages of < 



3 Myr due to the lack of evolution in colors a nd ionizing fluxes. 

an extinction value of Ay = 7.6 for WS80, and IGilbert & Graham! 1 

Ay = 10 for a simple screen mo del. Interestingly, in order to account for the lack of detection of 



Whi tmore & Zhang! d2002h infer 



(120071) find A K = 1 .2, or roughly 



H 2 0-SE at shorter wavelengths. ISnijders et al.l (|2006l) suggest that A v > 72 for H 2 0-SE! Given its 
very high extinction, H 2 0-SE may be an extremely young super star cluster in the making, though 
it is difficult to distinguish between extreme youth (and thus small physical size and Q LyC compared 
to WS80) vs. a simply less massive (and thus less able to clear its natal material) but more evolved 
cluster. 
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4.2. H 2 0-West Region 

The western maser resides in an optically thick dust lane coincident with a CO(3-2) molecular 
cloud, and is virtually invisible at every other available wavelength except for weak radio emission 
(see Fig. |2})). Moreover, because the radio emission is dominated by non-thermal emission it is 
impossible to discern the nature of any ionized thermal emission that might be associated with the 
maser region. It is possible that this source would be apparent at mid- to far-IR wavelengths, but 
unfortunately none of the available data have sufficient resolution to identify it. It is also notable 
that this maser has an exceptionally large luminosity (7.7 L ) for a GA-H 2 maser. Hopefully 
future high resolution near to mid-IR data may help shed light on this deeply embedded maser 
source. 



5. Conclusions 

We have imaged the GBT-discovered GA-H 2 masers in the Antennae galaxies with exquisite 
angular resolution (~ 80 mas), and find two maser regions in the IAR and a third 5"6 (> 600 pc) 
west of the NGC4039 nucleus. All three masers show excellent kinematic and spatial agreement 
with dense CO molecular gas. The IAR maser regions are located in areas of high optical/near-IR 
extinction, and are coincident with thermal ionized gas suggesting the presence of several thousand 
O*. Indeed, both of these maser regions seem to pinpoint extremely young sites of deeply embed- 
ded super star cluster formation. These results highlight the promise of using GA-H 2 masers 
to precisely locate the earliest phases of extragalactic cluster formation. The current isotropic lu- 
minosity sensitivity limit of 0.6 L , is still brighter than most Galactic H 2 masers, suggesting 
we may only be seeing the bright tip of the maser distribution. With its new broad bandwidth 
correlator, the Expanded Very Large Array will be an ideal instrument for future GA-H 2 maser 
studies (the 1.3 cm data presented here could be taken in 1/3 the time, while achieving 10 x better 
continuum sensitivity). 

This research used archival Spitzer and Hubble Space Telescope data, operated by the Jet 
Propulsion Laboratory and Space Telescope Science Institute, respectively under NASA contracts. 
We thank B. Brandl, C. Wilson, D. Iono, and G. Petitpas for providing their published data in 
digital format. K. E. J. acknowledges support from NSF through CAREER award 0548103 and 
the David and Lucile Packard Foundation through a Packard Fellowship. 
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